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We show that the recently reported anomalies in 6 —» sg + g~ transitions, as well as the long¬ 
standing — 2 discrepancy, can be addressed simultaneously by a new massive abelian gauge 
boson with loop-induced coupling to muons. Such a scenario typically leads to a stable dark matter 
candidate with a thermal relic density close to the observed value. Dark matter in our model 
couples dominantly to leptons, hence signals in direct detection experiments lie well below the 
current sensitivity. The LHC, in combination with indirect detection searches, can test this scenario 
through distinctive signatures with muon pairs and missing energy. 


I. INTRODUCTION 

The lack of an acceptable Dark Matter (DM) candi¬ 
date within the Standard Model (SM) is a pressing phe¬ 
nomenological motivation for the existence of new physics 
(NP) beyond the SM. Dark Matter may very well be close 
to the weak scale, emerging from theories addressing the 
electroweak (EW) hierarchy problem of the SM. How¬ 
ever, these theories generally predict a variety of new 
phenomena and particles around the TeV scale, whose 
existence still remains to be established experimentally. 
Despite the lack of a NP discovery, a few measurements 
are in mild tensions with SM predictions. In particular, 
there is a growing array of anomalies involving muons 
IUEZ1- Moreover, DM searches at direct detection exper¬ 
iments strongly limit its interactions with quarks (see 
e.g. REF. 0), suggesting that DM around the EW scale 
might preferentially couple to leptons. In this paper, we 
thus take a data-driven approach, entertaining the possi¬ 
bility that the observed muon-related anomalies are the 
first signals of leptophilic DM at the EW scale. 

The list of anomalies involving muons starts with the 
long-standing puzzle of the anomalous magnetic moment 
of the muon, = (g^ — 2)/2. The BNL measurement |T 
exceeds the SM prediction by about 3 standard devia¬ 
tions [8], 

Aa M = a® xp - a® M = (287 ± 80) x 10" 11 . (1) 

More recently, the LHCb experiment reported on a series 
of anomalies in (semi-) leptonic B meson decays, which 
together point to a possible new source of b —>■ sg + g~ 
transitions at short distances. The perhaps most tan¬ 
talizing deviation occurs in the ratio Rk = B(B + —► 
K + g + g~)/B(B + —► A' + e + e“), which was observed 
about 2.6 a below the theoretically clean SM prediction, 
Rf 1 - 1 ~ 10~ 4 0, 

Rk p = 0.745±^(stat) ± 0.036(syst). (2) 

Furthermore, the measured decay rate for B° — ► 
K°*g + g~ was found to exceed the SM prediction in a 
particular region of phase space [2], a result supported 


by the latest LHCb data 0 . The observed decay rate for 
B® -A g + g~ is slightly below, but compatible with the 
SM prediction [7]. While none of these LHCb anomalies 
by themselves are significant enough to claim a discov¬ 
ery, it is intriguing that they point to a common Lorentz 
structure when interpreted as a signal of NP [SHUj • A si¬ 
multaneous explanation of A a^, however, requires more 
sophisticated model building. In this paper, we propose 
a simple toy model, which addresses both the — 2 and 
the various LHCb anomalies. 

Most NP interpretations of the b -A sg + g~ anoma¬ 
lies postulate the existence of a new state in the range 
of A ~ 1 — 10 TeV with tree-level couplings to muons 
and quarks, for instance a Z' gauge boson ® mm 
or a scalar leptoquark HMD- The same interactions 
typically contribute to — 2 at the one-loop level, but 
yield too small a contribution to explain the discrepancy 
in Eq. (TO), due to the suppression by the high scale, 
oc mzj. A 2 [5] . A way out is to generate the coupling of 
Z' to muon pairs only radiatively , so that contributions 
to b — > sg + g~ transitions and — 2 are both induced at 
the one-loop level by NP around the EW scale. As we ar¬ 
gue in this work, this requires a richer NP sector with an 
electrically neutral state, which is stable if the tree-level 
Z'g + g~ coupling is forbidden by a (spontaneously bro¬ 
ken) symmetry. Hence, addressing the aforementioned 
muon-related anomalies generally yields a DM candidate, 
which, by construction, is mostly leptophilic. The same 
NP interactions dominate DM annihilation in the early 
universe. We will demonstrate that a minimal model 
with the above properties typically leads to a stable DM 
candidate with a thermal relic density of the order of the 
observed value. 

The remainder of the paper is organized as follows. We 
introduce our phenomenological model in SEC. [Tl] and 
identify the parameter space to simultaneously accom¬ 
modate the g^ — 2 and b -A sg + g~ anomalies in SEC. |III| 
We then discuss the implications of collider constraints in 
SEC. |IV| T he resulting DM phenomenology is analysed 
in SEC.|V|and SEC.|Vl] We conclude and give an outlook 
on future experimental tests of our model in SEC. |VII| 
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II. A MINIMAL MODEL 

We consider the extension of the SM by a new (dark) 
sector, consisting of heavy leptons L , L c and quarks Q , 
Q c with vector-like gauge couplings, as well as two com¬ 
plex scalars, </> and y. Interactions with the SM are medi¬ 
ated by a new gauge boson Z' associated with an abelian 
U(l)x symmetry, under which only the particles of the 
dark sector are charged. The quantum numbers of all 
new particles are listed in TAB. [I] jSS], Since no chi¬ 
ral fermion carries U(1 ).y charge, the model is anomaly 
free [115] ■ Notice that our choice of U(l)x charges for¬ 
bids tree-level Z' couplings with SM leptons. This is 
crucial in order to simultaneously address the g^ — 2 and 
b — > sy , + anomalies with NP around the weak scale. 

Besides canonical kinetic terms, the relevant new in¬ 
teractions in the Lagrangian are 

£np D eB^ v X^ - X x h\x\ 2 \H \ 2 - \^ H \ct>\ 2 \H \ 2 

-V((/),x)- [y (lL)x + W (qQ)(j) + h.c.] , (3) 

where 

V(<f>,x)= ( r< / ) X 2 + h.c.) + A^|0| 4 + A x |y| 4 , (4) 

and l T = (vcl , q T = ( ul,cLl) T and H are the SM 
lepton, quark and Higgs doublets, while B ^ and X/ v are 
the hypercharge and U(l)x held strength tensors, respec¬ 
tively. The U(1 )a symmetry is spontaneously broken by 
the vacuum expectation value (vev) of </>, leading to a Z' 
mass of mz' = 2 g'{(/>)■ It furthermore lifts the mass de¬ 
generacy between the components of y = (y 0 + ix')/V% 

by 


5 = 


to : 


a _ 1 = _ 2 lM. 


mt 


(5) 


We further assume that the singlet y is inert, i.e. does 
not develop a vev. There thus remains an exact Z 2 sym¬ 
metry, under which only y and L are odd, which ensures 
that the lightest state of the spectrum is stable. For def¬ 
initeness, we choose r < 0 and mz, > m Xo , so that the 
scalar component of y, yo, is a DM candidate [55]. As 
we will argue in the following sections, the above frame¬ 
work accommodates the g^ — 2 and b —> sy + yT anomalies 
simultaneously, without conflicting with current collider 
constraints, only in particular regions of parameter space, 
where the DM candidate yo is around the weak scale and 
relatively lighter than the other dark sector states, 


TABLE I: New fields and their quantum numbers. All SM 
fields are neutral under U(1 )a". 



spin 

SU(3) C 

SU(2)i 

U(l)y U(l)x 

L, L c 

1/2 

1 

2 

-1/2 1 

Q, Q c 

1/2 

3 

2 

1/6 -2 

<t> 

0 

1 

1 

0 2 

X 

0 

1 

1 

0 -1 


ensures that the lightest dark-sector state is a spinless 
SM-singlet DM candidate. 

We now turn to the new interactions in Eq. ©• The 
first term describes kinetic mixing between the hyper¬ 
charge U(l)y and the dark U(l)x held tensors, which is 
strongly constrained by electroweak precision measure¬ 
ments J2D] and collider searches [2T]. The second term, 
the scalar Higgs portal, is bounded by direct detection 
and collider searches [22]. Since both interactions are 
not relevant to our discussion, we set e = A x h = 0 [50] ■ 

The last two terms mediate interactions between the 
dark sector and the SM. The Z' couples with strength g' 
to the current j p = + Sj p , where j q is the tree-level 

U(1).y current and 

Sf = £ PfAi (hYf' L ) + h.c. (7) 

U ' 1 + 5 "’ 

is the part of the current induced after U(1).y breaking. 
Here q 2 is the squared Z' momentum, f,f = £,vz,u,d, 
and 5ff = 1 for / = /' and zero otherwise. SM quarks 
mix with dark quarks Q after U(l)x breaking, while SM 
leptons and dark leptons L do not mix due to the different 
charge assignments under U(l)x (see TAB. [I]). The Z' 
coupling to SM quarks thus arises at tree level through 
the mixing w(qQ)((f >) + h.c.. It is (^-independent and 
isospin universal at leading order. In contrast, the Z 1 
coupling to leptons is radiatively induced at the one-loop 
level through the exchange of dark-sector fields via the 
interaction y(lL) y+h.c., see FIG.fH In the limit <C 
m x,Li which we envisage here, the lepton form factors 
satisfy T«(g 2 ) ~ T„ eI/e (q 2 ), since the isospin components 
of L = ( L°,L ~) are mass-degenerate at leading order. 
Calculating the one-loop diagrams of FIG. [l] in the limit 
q 2 = 0 yields 


r«(o) 


\y \ 2 

32t r 2 


Fz’{r, S ), 


( 8 ) 


m xo ~ m x‘ ~ m Q ~ m z ' ; (6) 

while the dark lepton mass tol may be either close to 
the DM state or well above it. We do not address the 
dynamical origin of such a spectrum, but simply achieve 
it by tuning the bare mass squared of y and r{</>). We 
use 1/5 as a measure of this tuning, which remains mild 
in the parameter region of interest. Besides addressing 
the anomalies, this small hierarchy of the spectrum also 


where r = m 2 Jm 2 n and the loop function Fz> is given 
in Eq. (A2). 

Notice that Fz>{t,5) vanishes in the limit of unbroken 
U(1 )jy, 5 —► 0, where scalar and pseudo-scalar compo¬ 
nents are degenerate. This stems from the fact that the 
form factor F«(0) formally arises from the local operator 
(0* D p 4>)(Il1 p £l) + h.c. after U(l)x breaking, induced 
by the one-loop diagram in FIG. [l] with two insertions of 
the vev {(/). 
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FIG. 1: Dominant one-loop contribution to the operator 
(<l>*D p c/))(p,L'y p Hl) yielding F )J/i (0) ^ 0- The Z' line is under¬ 
stood to be attached wherever possible. A similar diagram 
with /r —t Up and L~ —¥ L° induces the Z’ coupling to neu¬ 
trinos. 


As of the flavor structure of U(l)x interactions, we 
follow a phenomenological approach and introduce only 
interactions that are required to explain the observed 
anomalies. We thus assume that the dark sector only cou¬ 
ples to the left-handed fermions ml> and &lSl + sl&l. 
The magnitude of the Z'b^SL coupling is experimentally 
constrained from B s meson mixing. Allowing for a NP 
contribution of 0(10%) to the mass difference A Mb, 
leads to the bound 0 [33 6Tj 

|r ‘-l s 2 - 4 * 10 " (sKsv) ■ < 9 > 


FIG. 2: One-loop NP contribution to — 2 and the Z cou¬ 
pling to muons. The photon 7 is to be attached to L~, and 
the Z boson to either of the fermion lines. 



FIG. 3: Leading NP contribution to the Wilson coefficients 
Cg 10 . The shaded disk denotes the radiatively induced Z' 
coupling to muon pairs at zero momentum exchange, as shown 
in FIG. □ 


We do not address here the origin of these peculiar flavor 
structures in the lepton and quark sectors. We acknowl¬ 
edge, however, that explicit flavor completions of this 
model could lead to correlated effects in other meson- 
physics and leptonic observables [HIISHSZ]- 


III. EXPLAINING MUON-RELATED 
ANOMALIES 

New contributions to the muon anomalous magnetic 
moment are induced at the one-loop level through the 
diagram in FIG. [2j yielding 


the effective Hamiltonian 

-Hel = - ^ V tb V* OiOi + h.c., (11) 

v i 

where a and Gf are, respectively, the fine-structure and 
Fermi constants, Vij are CKM matrix elements, and the 
sum runs over the operators (£ = e, m) 

O e 9 = & 7 P (1 — 7 5 )s £"f p £, ( 12 ) 

C’io = h P (X -7s)s h p ld- (13) 

A global fit to leptonic, semi-leptonic and radiative B 
decays favors effective couplings to muons 116, 28] 


a < p = !&£fu T) i i+A(T - 4) i’ (io) 

with nip the muon mass and A(t,S) = Fp ifI)f + (t) • The 
loop function F g is given in Eq. ( |Al[ ). For a degenerate 
dark spectrum, we have F g (l) = 1/12 and A(1,0) = 1, 
so that the discrepancy Ain Eq. 0 is accommodated 
for m Xo ~ 451 y GeV. As we will see, the LHCb anomalies 
typically require a significant scalar versus pseudo-scalar 
mass splitting i > 1, with a dark-lepton versus DM mass 
splitting roughly in the interval r £ [1,5]. Eq. (TTl) is then 
accommodated for m Xo ~ 32|y|(l + 2/5) GeV (r ~ 1) or 
m Xo ~ 55|y|/V5GeV (r ~ 5 » 1). (See APP. |a| for 
details.) 

Consider now the b —t sg + ^~ anomalies. At the 
b— quark mass scale, the NP amplitude is described by 


eg = -CTo ^ -0.5 (14) 


and negligible electron coefficients Cf 10 ~ 0 [62]. In our 
model, C 910 are induced at the one-loop level through 
the diagram in FIG. [3] which gives 


9 _ ° 10 _ 4 ml 


9' 2 a 2 WL r 

iV SM Tr T 7" * 1 fcs 1 


VtbV* 


MM > 


(15) 


where Asm = [2\/2n/(aGF\V t bV t * I)] 1 / 2 — 50TeV is the 
scale of the SM contributions. Accommodating the b —► 
SM + M _ anomalies as in Eq. |l4| ), while respecting the 
A Mb s bound from Eq. yields a lower bound on the 
muon form factor, 


r w K) > 0.029 


mz'/g' 


300 GeV 


(16) 
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Since the Z' coupling to muons is small, cf. Eq. ([8]), this 
bound is fulfilled only for a light Z' boson, as well as 
a large mass splitting S. In a spectrum with moderate 
tuning, S < 10, this leads to 


m z > < 110 — 270 GeV (^j ^ 

where the two mass values correspond to r — 8 and r ~ 
1, respectively. For mz> around the weak scale, LHCb 
anomalies thus require rather large couplings y and g', 
while a perturbative upper bound of mz' < 8.5 TeV (r — 
<5) or mz> < 20 TeV (r ~ 1) applies for g' = y = 47t. We 
discuss collider constraints on such a Z' in SEC. IIV1 To 
summarize, assuming perturbative couplings g' ~ 3 and a 
mild tuning of the dark spectrum, 8 — 10, both the g^ — 2 
and LHCb anomalies are accommodated for parameters 
values interpolating between the two limiting cases 

1) r ~ <5, \y\ ~ 6, g' ~ 3, 5 ~ 10, 

m Xo ~ 100 GeV, ~ 300 GeV, (18) 

2) r~l, |y| — 2, </~3, <5-10, 

m Xo ~ 70 GeV, m Z ' ^ 150 GeV. (19) 

IV. COLLIDER CONSTRAINTS 

We now analyse the relevant constraints on our model 
from EW precision measurements at the LEP experi¬ 
ments, as well as from the first LHC run. First of all, 
our model implies sizeable radiative corrections to the 
Zgp,, Zv^v^ and W + fj,u fl couplings from the one-loop 
diagram in FIG. [2] The EW gauge couplings are shifted 
by {V = W,Z) 



^-^Fyirr) 
- " 32^ 2 V[ ,rq) 


9sm 


q) > 


( 20 ) 


where the one-loop function Fy is found in Eq. (A4) 


and r q = q 2 /mr L . The vertex correction is my/mr L - 
suppressed at the V pole, q 2 = rriy. While the QED part 
of the couplings at zero momentum is protected by gauge 
invariance, the weak isospin part is not corrected at one- 
loop level, since the SM Higgs doublet does not directly 
couple to the dark sector JBB] - Despite this paramet¬ 
ric suppression, the lightness of the dark-sector states, 
m Xo — 100 GeV and tol — 100 — 400 GeV, together with 
a relatively large Yukawa coupling, y > 2, typically shifts 
the SM gauge couplings by one to a few permil, which 
is in mild tension (~ 1 <j — 3<r) with LEP data [B]. This 
tension may be relieved in a more sophisticated version 
of the minimal model considered here. 

Our model also predicts a series of signatures at hadron 
colliders, most notably muon pair production through a 
resonant Z', as well as signals of large missing energy 
with muon pairs and/or jets. Current LHC limits on 
a Z' resonance with SM-like couplings to fermions are 


around mz> > 3 TeV ESEHl- However, Z' production 
in our model only occurs through sea-quark (bs + sb) 
annihilations and is thus strongly suppressed. On the 
other hand, the Z' dominantly decays into muon pairs 
and neutrinos. Given the conditions Eqs. ([9]) and © 
on the Z' couplings to SM fermions and for Z' masses 
that accommodate the LHCb anomalies as in Eq. ( fl7| , 
the branching ratios to leptons are both ~ 40 — 50%, de¬ 
pending on the value of the Yukawa coupling (BJ). We 
thus find the cross section for Z 1 production with decay 
into to be of O(fb), which is an order of magni¬ 

tude below current limits at the 8 TeV LHC. The next-to- 
leading Z' branching ratio is g + g~XoXo, ranging from 2 
to 10% for large Yukawa couplings. Along the same lines, 
mono-jet signatures from the direct production of a DM 
pair in association with a hard jet from initial state ra¬ 
diation (ISR) lie at least one order of magnitude below 
the current LHC sensitivity EUEg. 

The EW production of L + L~ pairs leads to a signa¬ 
ture with di-muons and missing energy, which resembles 
the one used in searches for smuons, the supersymmetric 
partners of the muon. The only difference with our signal 
lies in the spin of the produced particles. However, it was 
shown in REF. [33] that results for slepton searches in 
simplified models could safely be applied to the produc¬ 
tion of fermion pairs decaying into a fermion and a scalar. 
We use SmodelS [33], a tool designed to decompose the 
signal of any NP model into simplified topologies, and 
compare the predictions to the exclusion limits set by 
the ATLAS and CMS slepton searches 35 35]. We find 
that the 8 TeV LHC sets strong constraints on the mass 
of L ~, even stronger than for smuons, because L + L~ 
pair production cross sections are significantly larger. 
Dark-lepton masses rn^ < 450 GeV are excluded, except 
if the mass splitting with the DM is sufficiently small, 
tol — m Xo < 60 GeV. In this region, the di-muon signal 
is overwhelmed with SM background. Similar searches 
at LEP2 lead to the lower bound of Tol > 100 GeV [ 37] , 


V. DARK MATTER RELIC ABUNDANCE 

The DM candidate in our model is the lightest com¬ 
ponent of the scalar Xi which we assume to be Xo- It 
is largely leptophilic, as follows from the charge assign¬ 
ments in TAB. [I] For a spectrum as in Eq. (|6|, DM 
annihilation proceeds dominantly into g + g~ and 
through t —channel exchange of L~ and L°, respectively, 
as shown in FIG. [4j Coannihilation processes with the 
other dark states are negligible for r ~ S 1, while anni¬ 
hilation into Z' pairs is negligible as long as mz> > m Xo . 
The resulting annihilation cross section is d —wave sup¬ 
pressed in the chiral limit, ayv ~ v 4 [38] . Adding up 
final-state muons and neutrinos, the thermal average 

is [35] 


{°Tl v ) 


ad 

x 2 ml 0 


+ 0(x ~ 3 ), 


„ _ \y\ 4 

d 27t(1 + t) 4 ’ 


( 21 ) 
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FIG. 4: Dominant amplitudes for DM-DM annihilation (top) 
and DM-I or L-L coannihilation (bottom) in the early uni¬ 
verse. 


where v is the relative DM-DM velocity, and x = m Xo /T, 
the DM mass-to-temperature ratio. For the parameter 
set in Eq. (181, ad — 0.02 is a typical value. In the 

, the relic density is 


freeze-out approximation 
-i-ii I 


H“/i 2 ~ 8.5 x 10" 


3a V w x 0 


Ud^GeV 2 


x [l + R f }-\(22) 


where g * ~ 80 is the effective number of relativistic de¬ 
grees of freedom at T/ ~ 1 — 10 GeV and 


R f ~ 2.6 x 10 


— 18 ( V9* m Xo \ 3/2 


V a d GeV J f 


x / exp(xf) (23) 


accounts for the reduced DM abundance at freeze-out. 
Since DM annihilation (agv) oc x ~ 2 is d— wave sup¬ 
pressed during freeze-out, DM decoupling is slightly de¬ 
layed, which results in a larger (power of) Xf ~ 25 in 
Eq. (23) compared to s— or p —wave annihilation. In a 


scenario as in Eq. (18), Rf ~ 0(1) is a significant cor¬ 
rection. 

Without this correction and for fixed mass splittings 
r and 5, the relic density scales as f l l x h 2 oc m Xo /\y\ i , 
i.e. like (Ax Cg) -1 , in our model. Using Eqs. (fTol) 
and (|l5|, Eq. (22) reads (Rf = 0) 


'x f \ 3 /100 GeV\ 


Sl“ft 2 =2 0.01x/(r,i)(|i)‘ 


287 x 10 
A a„ 


-li 


£)• <“) 


where 1(6, t) = F g (r) [1 + A(t, <5)] F z >(t, <5)(1 + r) 4 . As¬ 
suming parameter values of Eq. ( fl8] ) (with r 1) 
in order to accommodate the muon anomalies yields 
J(r, <5) ~ 20 — 40. Recalling that Rf ~ 0(1), this gives a 
relic abundance close to the observed value from Planck 
data H01 


D° bs /i 2 = 0.1199 ±0.0022. 


(25) 


However, for parameter values of Eq. (19) (with r ~ 1), 


J(r, 6) ~ 0(1) and the resulting relic density is typically a 
factor of 5 —10 smaller than observed. We stress that the 
d —wave suppression of scalar DM annihilation into lep¬ 
ton pairs through a heavy fermion mediator plays an im¬ 
portant role in the above prediction. Had DM annihila¬ 
tion been p —wave (s— wave) dominated, the same param¬ 
eter space would have predicted a relic density smaller by 
a factor of about ten (hundred). 


Since the two-body process in Eq. (21) is strongly 
suppressed [66] . parametrically sub-dominant annihila¬ 
tion channels can become relevant. The leading such 
contribution is three-body annihilation XoXo -± A i+ M _ 7 
with an extra photon in the final state. Virtual inter¬ 
nal Bremsstrahlung (VIB) from the dark lepton L + lifts 
the kinematic suppression, so that annihilation proceeds 
through an s —wave [38, 41]. The thermally averaged 
cross section for this process is [42] 


yV) = 


\y \ 4 


327r 2 ' 


■F 7 (t) 


(26) 


Xo 


with the function F 7 defined in Eq. (|A6|. Following 
REF. [25], we find that VIB suppresses the predicted relic 
density in Eq. (22) by a factor i? 7 = < 1, 


which in the chiral limit and at zero relative DM velocity 
v = 0 approximates m 

~ 1 ++ ■ ( 27 ) 

The function (1 ± t) 4 F 7 (t) strongly decreases with r, 
yielding to suppression factor of I? 7 ~ 0.5 (0.9) for r ~ 1 
(r> 1). Therefore, VIB is a significant contribution only 
in particular regions of parameter space, where coanni¬ 
hilation in addition strongly depletes the relic density. 

In order to assess the robustness of the prediction for 


the relic density in Eq. (24) more quantitatively, we 


broadly scan around the parameter sets in Eqs. (18) 


and (19) that accommodate the muon anomalies, varying 


the parameters in the following ranges 


50 < m Xo < 450 GeV, 200 < m z > < 1000 GeV, 

1 < r < 20, 0< 5 < 20, 1 < |y| < 4 tt, 1 < g <5.(28) 

We restrict the tuning arnoung the scalar components to 
maximally 10%. Furthermore, we impose the condition 
that the g^ — 2 and the LHCb anomalies quoted in SEC. [T] 
are explained within one standard deviation, and that the 
collider limits on di-muon and missing energy produc¬ 
tion from SEC. |IV| are satisfied. We then compute the 
resulting relic density with micrOMEGAs [23]. As shown 
in FIG. [5] for many parameter points where the domi¬ 
nant DM annihilation process is XoXo U (dark gray 
points) the relic density falls within an order of mag¬ 
nitude of H° bs , as expected from the above discussion. 
However, we also find that the predicted relic density 
can be much lower than the observed value, if the spec¬ 
trum significantly deviates from Eq. § First of all, 
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f l x rapidly decreases when the DM mass increases above 
ro Xo > 200 GeV. Satisfying Awith heavier DM masses 
requires a lighter dark lepton L or, equivalently, smaller 
values of r, leading to a significantly stronger DM annihi¬ 
lation into muon(-neutrino) pairs, see Eq. (21). Secondly, 
when mz' < m Xo , the process XoXo —1 > ZTZ' efficiently 
depletes the DM relic abundance (green points in FIG. [5]) . 
Moreover, when rriL — rn Xo , corresponding to r ~ 1, 
coannihilation processes with dark leptons (see FIG. [4]) 
can also contribute to significantly reducing the DM den¬ 
sity (magenta points in FIG. [5]). This is particularly 
pronounced when Z' can be produced in the final state 
via the coannihilation processes Xo ^(L 0 ) —>• Z' ^(y^) 
(cyan points in FIG. [5]). 

The two regions where the relic density is in agree¬ 
ment with the observed value closely resemble the param¬ 
eter limits favored by the muon anomalies from Eqs. (18) 
and (191, respectively. These are 


1) The heavy lepton scenario. In this region, m Xo ss 
120 — 250 GeV and mi > 450 GeV, corresponding 
to r > 1, because of collider constraints on di¬ 
lepton plus missing energy signals. As a result, the 
Yukawa coupling is strong, \y\ > 6, to ensure a large 
enough DM annihilation rate into lepton pairs. 


2) The compressed scenario. This second region cor¬ 
responds to m Xo « 50 — 80 GeV and tul ~ 100 — 
120 GeV, i.e. to r ~ 1. It evades collider con¬ 
straints because of the relatively small gap between 
DM and dark-lepton masses. Since DM annihila¬ 
tion is larger for r ~ 1, see Eq. (21), a smaller 
Yukawa coupling is needed to ensure enough DM 
annihilation. However, as the muon anomalies re¬ 
quest | y | > 2, it is more challenging to obtain the 
observed relic density 0“ bs h 2 ~ 0.1 in the com¬ 
pressed scenario. 


In FIG. [6j we illustrate the range of y needed to obtain 
various amounts of DM relic abundance f2 x , as a function 
of the DM mass, for the case of dominant DM annihila¬ 
tion into lepton pairs, XoXo —> U- 


VI. DARK MATTER DETECTION 

The DM candidate in our model is dominantly lep- 
tophilic, which makes direct detection very challenging. 
Furthermore, since Z and Higgs boson couplings to xo 
pairs are absent at tree level and loop-induced vector in¬ 
teractions are momentum-suppressed, any signal of spin- 
independent DM-nucleus scattering lies below the sensi¬ 
tivity of current direct detection experiments [68]. 

On the other hand, VIB in DM annihilation, discussed 
in SEC.[V] plays an important role for indirect detection 
today. The suppression of tree-level DM annihilation into 
di-lepton states, ayv ~ v 4 , is even stronger than during 
freeze-out, as velocities today are around I'haio ~ 10 -3 in 
our galactic neighborhood. Photons from VIB are thus 



m x o [GeV] 


FIG. 5: DM relic density Q l x h 2 as a function of the DM mass, 
m xo , as predicted by the muon-related collider anomalies. 
Light gray points: excluded by 8TeV LHC and LEP2 con¬ 
straints on di-lepton plus missing energy signals from L + L~ 
production Dark gray points: DM annihilation dominated by 
XoXo —> p + pTfVpVy.. Green points: additional contributions 
(~ 10% and more) from XoXo —> Z'Z' . Magenta and cyan 
points: co-annihilation through LL —> ll and xo L —> Z'l, re¬ 
spectively, accounts for more than ~ 10% of the total (crv). 
Red band: 3<r range of S2 x bs from Planck [40]. (See text for 
details.) 



100 200 300 400 

m Xo [GeV] 


FIG. 6: Dark Yukawa coupling y as a function of m xo re¬ 
quired to accommodate the muon anomalies within one stan¬ 
dard deviation. Gray points are excluded by 8TeV LHC and 
LEP2 constraints on di-lepton plus missing energy signals 
from L + L~ production. All other points correspond to dom¬ 
inant DM annihilation XoXo H ► y + y~, i.e. to the dark 

gray points in FIG. [5] Colors depict different amounts of the 
associated DM relic abundance; Q x h 2 is obtained in the blue 
regions. 


the dominant signal in gamma-ray searches of indirect 
detection experiments. The annihilation cross sections 
due to VIB in our two scenarios from Eqs. (18) and (19), 
integrated over the photon energy spectrum, are 


1) ~ 3 x 10 27 cm 3 /s and 

2) ( W ) ~7xl0- 25 cm 3 /s, (29) 


respectively. The smaller cross section in the heavy 
lepton scenario 1) is mainly due to the suppression 
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FIG. 7: Thermally averaged cross section for DM annihilation 
into fi + fj,~'y in our galactic neighborhood today as a func¬ 
tion of the DM mass, as predicted by muon-related collider 
anomalies. All points satisfy 0.01 < ffy/i 2 < 1. Light gray 
points are excluded by 8TeV LHC and LEP2 constraints on 
di-lepton plus missing energy signals from L + L~ production. 
Points with r < 2 (r > 2) are shown in blue (dark gray). The 
black line denotes the Fermi LAT upper bound on photon 
emission in nearby dwarf galaxies, while the dashed red line 
represents the Fermi LAT limit on spectral features in VIB 
signals j3U, which applies only for 1 < r < 2, i.e., to the blue 
points. 


( cr fip,'i v ) ~ l/ 7-4 hi the limit of large mass splitting, r 1. 

For small mass splitting r ~ 1, the energy spec¬ 
trum of VIB photons is peaked close to the endpoint 
Ery = m Xo |41j . which leaves a characteristic signature 
in indirect detection experiments. The Fermi satellite 
has performed a dedicated search for spectral features 
in VIB using data from space regions near the galactic 
center collected during Pass 7 [44] . Assuming a spectral 
shape characteristic for 1 < r < 2, this search sets an 
experimental upper bound on VIB, 

(cr^v) < 10 _27 cm 3 /s, (30) 


for m Xo < 100 GeV. 

More inclusive analyses of Fermi data from dwarf 
galaxies constrain the overall gamma-ray flux emitted in 
DM annihilations [IS] [3B] . These searches are sensitive 
to the sum of secondary photons (from final-state radia¬ 
tion off muons and muon decays) and VIB. Though the 
cross section bounds are generally weaker than from the 
dedicated VIB search, they constrain (cr M ^ 7 u) for scenar¬ 
ios with large r, which the VIB spectral search is not 
sensitive to. 

In FIG. [T] we compare the predictions of DM annihi¬ 
lation through VIB in our model with the bounds from 
Fermi. We focus only on points in the scan range of 
Eq. (281, which satisfy the muon-related anomalies and 
yield a DM relic density of 0.01 < f l x h 2 < 1. The upper 
bound on (<T fl p, 1 v) from Fermi’s VIB search was derived in 
REF. [44] (dashed red line). It excludes part of the com¬ 
pressed scenario 2) with 1 < r < 2 (blue points), where 
VIB is particularly large (see Eq. (29)). The bound can¬ 
not directly be applied to scenarios with r > 2 (dark 


gray points), for which the photon energy peak is much 
less pronounced and the VIB search thus significantly 
loses sensitivity. Inclusive gamma-ray searches in nearby 
dwarf galaxies lead to a less strigent upper bound on 
(black line), which we adopt from [37] under the 
assumption that VIB dominates photon emission in DM 
annihilation. Parameter space regions with r > 2 and 
a relic density comparable to the observed value lie at 
least an order of magnitude below the bound from dwarf 
galaxies. Notice that the Fermi bounds on the cross sec¬ 
tion are based on data from Pass 7. With the recently 
released Pass 8 data set [3B], we expect these bounds to 
improve by about an order of magnitude. The inclusive 
gamma-ray search might thus start probing the heavy 
lepton scenario, as well as the part of the compressed 
scenario, which does not display a significant spectral 
feature. 


VII. CONCLUSIONS AND OUTLOOK 

We have shown that a minimal phenomenological 
model featuring a new abelian force with radiative cou¬ 
plings to left-handed muons can simultaneously explain 
the observed muon anomalies in — 2 and b —» sp + 
transitions through one-loop effects of new particles 
charged under the associated U(l)x symmetry. The 
spontaneous breaking of U(l)x down to a Z 2 symme¬ 
try induces a sizeable mass splitting among the scalar 
states of the dark sector, which controls the radiative Z' 
coupling to muon pairs and ensures that the lightest new 
state is a cosmologically stable neutral scalar. Hence, 
addressing both the — 2 and b —» s[x + anomalies 
in this way typically yields a leptophilic DM candidate, 
which is hardly visible in direct DM detection searches. 

All muon-related anomalies can simultaneously be ac¬ 
commodated with Z' and DM masses around the EW 
scale and large dark couplings y, g' > 2. While DM anni¬ 
hilation is governed by the same interactions that explain 
the muon anomalies, the large Yukawa coupling y typi¬ 
cally allows for a thermal DM relic density surprisingly 
close to the observed value, unless significant coannihila¬ 
tion is present. This result partially relies on the fact that 
the dominant DM annihilation into lepton pairs is d -wave 
velocity-suppressed, thus favoring larger couplings com¬ 
pared to scenarios with s— or p— wave annihilation. The 
velocity suppression furthermore depletes final-state and 
secondary photon emission from DM annihilation into 
muon pairs in our galactic neighorhood today. Signals in 
indirect detection experiments are therefore dominantly 
due to hard photons from VIB. 

Another key ingredient of the model is a new vector- 
like lepton doublet, which in particular mediates tree- 
level DM annihilation to leptons. Being Z 2 -odd, these 
dark leptons do not mix with SM leptons, nor do they 
couple to the SM Higgs field at tree level. This prevents 
the otherwise severe constraints from Higgs and EW pre¬ 
cision observables, most notably oblique corrections. The 





only significant effect on low-energy observables lies in 
the Z boson coupling to muon pairs or invisible states, 
which is in mild tension with precision LEP data at the 
Z pole. 

At high-energy hadron colliders, dark leptons around 
the EW scale are produced in pairs through Drell-Yan- 
like processes. We recast available LHC searches for 
muon pairs with missing energy in the context of super- 
symmetry and find a general lower bound on the dark 
lepton mass, > 450 GeV. A good fit to all muon 
anomalies with a relic density around the observed value 
is still possible in two characteristic parameter regions. 
These are 1) the heavy lepton scenario with dark lep¬ 
tons above 450 GeV and a DM candidate in the range 
125 < m Xo < 250 GeV, and 2) the compressed sce¬ 
nario with a light dark sector around 70 GeV and a 
small splitting between the dark-lepton and DM masses, 
?7i£ — m xo < 60 GeV, to which current di-muon searches 
are insensitive. 

Accommodating the LHCb anomalies in the heavy lep¬ 
ton scenario typically requires a very large Yukawa cou¬ 
pling not far from the non-perturbative regime, which 
suggests a UV-completion of our model not far above the 
TeV scale. Weaker couplings are nevertheless possible in 
the compressed scenario thanks to the larger DM anni¬ 
hilation cross section at freeze-out in the limit of a small 
mass gap with the dark lepton. However, this small mass 
splitting induces a sharply peaked photon energy spec¬ 
trum from VIB in DM annhihilation today. Dedicated 
searches for such a signature by the Fermi experiment 
partially exclude the compressed scenario as a simultane¬ 
ous explanation of muon anomalies and thermal DM. In¬ 
terestingly, indirect detection experiments are thus com¬ 
plementary to collider searches, which are not sensitive 
to the compressed scenario. Probing the heavy lepton 
scenario through indirect detection, however, is difficult, 
since the photon energy peak from VIB is much less pro¬ 
nounced for large mass gaps and inclusive gamma-ray 
searches are less sensitive. 

We point out that B s meson mixing significantly lim¬ 
its the Z'hs coupling to at most a few perrnil. Given the 
loop suppression of the Z' coupling to muon pairs in our 
model, the b —> sg + anomalies imply new contribu¬ 
tions to B s meson mixing right around the corner, while 
they are only a possibility in most Z' models explaining 
the LHCb anomalies with tree-level couplings to leptons. 
This is an interesting implication of simultaneously ad¬ 
dressing the g g — 2 discrepancy. Moreover, the absence 
of a Z' coupling to valence quarks and its constrained in¬ 
teraction with bottom and strange quarks strongly limits 
direct Z' production at hadron colliders. In particular, 
resonant Z' production followed by decays in muon pairs 
is well below the sensitivity of the 8TeV LHC. Similar 
conclusions apply to mono-jet searches. 

Prospects are good that our model can be tested in 
the near future. Since dedicated DM searches through 
direct detection are not sensitive to our scenario, the 
14 TeV LHC will be the most important experiment 


to probe it. We expect that the reach for dark muons 
should be improved by close to a factor of two during 
run II, similarly to what was found for srnuons [151 US]- 
Thus, the parameter space of the heavy lepton scenario 
will be fully covered, since the — 2 anomaly sets an 
upper bound tol < 800 GeV for couplings y < 47r. The 
spectrum in the compressed scenario (which is already 
under pressure from indirect DM detection searches) is 
more difficult to access through direct searches. One 
possibility would be to exploit the monojet signature 
from the direct production of dark muons in association 
with an ISR jet in the highly compressed region, where 
the dark muon decay leads to an invisible signature. 
Precision measurements of Z boson couplings at a 
future e + e _ collider would give complementary indirect 
information about the dark sector. 
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Appendix A: Auxiliary Functions 


The g g — 2 loop function from the diagrams in FIG. [2] 
is 

F g( T ) = 6( - r 1 “ 6r ( r - lo S r ) + 3r + 2] ■ (Al) 

For degenerate dark (pseudo-)scalar and dark lepton 
masses, one has F g (r = 1) = 1/12, while in the decou¬ 
pling limit of the dark lepton the function decreases as 
F 9 {t » 1) ~ 1/(6t). 


For the dominant Z' coupling to SM lepton pairs 
at zero-momentum, corresponding to the diagrams in 
FIG. [T[ the loop function is 


F Z '(t,S ) 


1*1 1*1 —U ' 

1 Hi Hi 1 Hi H J- I / 

JO Jo L T+('U + ' 

1 Ion [l v6 1 

■ (A2) 

t + (w + u)(1 + <5-t)_ 


For degenerate scalar and pseudoscalar masses, 
Fz>(t, 0) = 0, while in the limit of large mass splitting 
Eq. (|A2|) approximately gives 


F Z ' (t, S » 1) 



r 2 ( 1 + logr) — 3r + 2 


.(A3) 


2( r — l) 2 
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Eq. (A2) approaches a constant when 
are comparably large, 


Moreover 
both mass splittings 
F Z /(r~<y» 1) ~ 1/4. 


The one-loop corrections to the V = W, Z couplings 
to SM leptons from FIG. [5] are proportional to the loop 
function 


f 1 

✓*1 — U 

/ du 

/ dv 

1 0 

Jo 


1 — u — v 
u + v 


log f(j,u + v) + 


1 — /(r, u + v) + 2 uvr q 
fir, u + v) — uvr q 


- log(/(r, u + v)- uvr q ) 


,(A4) 


where f(r,x) = x + (1 — x)/t and r q = q 2 3 /m 2 Ll with where Li 2 (a;) = — dy \og(l—xy)/y is the polylogarithm 

q 2 the squared momentum of the external gauge boson. function of second order. 


The loop function falls off at low q 2 as 

rqT [(r - l)(7r 2 - 29r + 16) 

(A5) 


Fv(T,r q < 1 ) ~ 


36(r - l) 4 5 * * 8 9 10 11 12 13 14 15 
+6(3r - 2) log r]. 


Finally, the function controlling the approximate 
3—body correction from VIB to the relic density is 


-F 7 (t) = (r + 1) 


~6 l0S “ 


r + 1 
2r 


— 2Lio 


4r + 3 4 r 2 — 3r — 1 


r + 1 


27- 


log 


r — 1 
r + 1 


r + 1 
2r 


, (A6) 
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